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Abstract. The axially symmetric vibrations of an isotropic poroelastic cylindrical panel 
resting on elastic foundations are investigated in the framework of Biot’s theory. The effects of 
the surrounding elastic medium are considered using the spring constant of the Wrinkler type 
and the shear constant of the Pasternak type. The frequency equation is obtained for both 
pervious and impervious surfaces. Non dimensional phase velocity is computed as a function of 
wavenumber. Numerical results are presented graphically. 


1. Introduction 

The study of elastic foundations in cylindrical structures has wide applications in the field of 
Engineering, Geology and Biomechanics. Cylindrical panel plays an important role in Aerospace 
as they are of high specific strength and high specific stiffness. Even in human body, a num- 
ber of muscoskeletal models of knee joint and skin employ different forms of elastic Wrinkler 
foundations. Propagation of axially symmetric waves in finite elastic cylinder is studied in [1]. 
In the said paper, vibrations corresponding to the different propagation constants for the same 
frequency parameter are applied to satisfy the boundary conditions. In [2], authors studied wave 
propagation in a homogenous isotropic cylindrical panel embedded on an elastic medium. In 
the said paper, three dimensional wave propagation in cylindrical panel embedded in an elastic 
medium (Wrinkler model) is investigated in the context of linear theory of elasticity. Free vi- 
brations of circular cylindrical shell on Wrinkler and Pasternak foundations is studied in [3]. In 
said paper, it is concluded that elastic foundation affects radial vibrations mode frequency while 
pertaining to torsional and longitudinal remains are unaffected. Employing Biot’s theory [4], 
axially symmetric vibrations of a poroelastic composite solid cylinder is investigated in [5]. In 
the said paper, frequency equation is obtained for nonaxially and axially symmetric vibrations 
each for pervious and impervious surfaces, and the results are compared with that of rule of 
mixture (RoM). In [6] authors investigated vibration characteristics of fluid filled cylindrical 
shells based on elastic foundations. In the said paper, frequencies are strongly affected when a 
cylindrical shell is attached with elastic foundation. Axially symmetric vibrations of fluid filled 
poroelastic circular shells is studied in [7]. In [8], authors investigated axially symmetric vibra- 
tions of finite composite poroelastic cylinders. In the said paper, non-dimensional phase velocity 
for propagating modes is computed as a function of ratio of length of cylinders in the absence 
of dissipation. To the best of knowledge, axially symmetric vibrations in poroelastic solid cylin- 
drical panel resting on elastic foundation are not yet investigated. Therefore, in this paper the 
same is investigated in the framework of Biot’s theory. Frequency equations are obtained for 
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both pervious and impervious surfaces. Non-dimensional phase velocity against the wavenumber 
for different elastic foundation are computed for three solids namely sandstone saturated with 
kerosene, water, and bone. 


This paper is organized as follows. In section 2, governing equations, and solution of the 
problem are given. In section 3, frequency equations are derived. Numerical results are done in 
section 4. Finally, conclusions is given in section 5. 


2. Governing equations and solution of the problem 
The equations of motion of a poroelastic solid [4] in presence of dissipation (b) which in terms 
of displacement vectors are 


NV? + (A+ N)Ve + Qe = & (put + pÜ) +088- Ü), 
(1) 


V(Qe + Re) = Z (pit + pÜ) — b2 (a - Ü), 


where, V? is the Laplace operator, (u, v, w) and Ü (U, V, W) are solid and fluid displacements, e 
and € are the dilatations of solid and fluid respectively; the symbols A, N, Q, R are all poroelastic 
constants; pij are mass coefficients. The constitutive relations are 


Ou = 2Neij + (Ae + Qe)dj;, (i,j = 1, 2, 3), 
(2) 
s = Qe + Re. 


In eq. (2), eij’s are strain displacements, o;;’s are solid stresses and fluid pressure s, 6;; the well 
known Kroneckar delta function. For axially symmetric vibrations, the displacements of solid 
u(u,0,w) and fluid U(U,0,W) which in terms of potential functions ¢’s and y's are given below 
by 


a ee ee res 
Ab, _ ay Ags , ðv, v K 
Uea = Bae Waa Fe eae 


For free harmonic waves travelling in the z- direction, we take 


pı = Fi(r)coskze™t, p2 = Fo(r)coskze™t, 


(4) 
yı = Gi(r)sinkze™, pz = Go(r)sinkze™. 


In eq.(4), k is the wavenumber. w is the frequency of the wave, i is the complex unity. By 
substituting eq. (4) in eq. (3), we obtain solid displacements as follows. 


u = —(CipKi (pr) + CoqKı (qr) + Aik Ki (dr))coskze™, 


(5) 
w = —(CikKo(pr) + CokKo(qr) + AidKo(dr))sinkze™*. 


In eq. (5) C1, C2, Ai are all arbitrary constants, Kn (x) is the modified Bessel functions of second 
kind of order n and p,q, d = k(1—€?), & = aoi = 1,2,3. Here V;(i = 1, 2,3) are the dilatational 
wave velocities of first and second kind and shear wave velocity, respectively. Making use of eq. 
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(2) and eq. (5), we obtain the relevant stresses and fluid pressure that are given below 
Orr +8 — Ku — GAu = (An (r)C1 + A12(r)C2 + A13(r) Ai)coskze™, 
orz = (Aai(r)C1 + A22(r)C2 + Ao3(r)A1)sinkze™, 
s = (Asi (r)C1 + Aga(r)C2)coskze™, 
gs = (B3i(r)C1 + B32(r)C2)coskze™, 


where, 


= 
Í 


2NP Ky (pr) + (2Np? + (P — 2N + Q8?) (p? — k?) + (Q + Rô?) (p? — k?) Ko(pr) — 7%? Ki (pr), 
Ais = 222 Kı (dr) + 2NkKo(dr) — Z Kı (dr), 

Agi = 2NkpKi (pr), 

A23 = N(d? = k?) Ki (dr), 

Agi = (Q + Rô) (p? — k?)Ko(pr), 

A33 = 0, 


oe (PR-Q?) VT ?-(Rm11-Qmı2) 
Le Rmi2—Qm22 ? 


(7) 
63=similar expression as ô? with V; replaced by V2,A12, 422, 432=similar expression as A11, A21, 
A31,with p and 6?,replaced by q and 63,respectively, and m11 = p11 — ibw—!, m12 = pig + ibw7!, 


_ aes 
M22 = p22 — ibw™ `. 


3. Boundary conditions and frequency equation 
The boundary condition for stress free surface at r = a for the pervious surface is 


Orr + s — Ku — GAu = op = 8s =Q. (8) 


The boundary condition for stress free surface at r = a for an impervious surface is 


Orr +5 — Ku — GAu = or = & = 0. (9) 
In eq. (8)and (9), A = & +42 + cs K is the foundation modulus given by K = Taya? where 


E is the Young’s modulus, p is the poisson’s ratio [9]. If the shear modulus G goes to zero, the 
Pasternak foundation will reduce to Wrinkler foundation [10]. Eq. (8)and (9) gives a system 
of three homogenous equations in the three arbitrary constants C1, C2, A; each for a pervious 
surface and impervious surface. A nontrivial solution can be obtained if the determinant of the 
coefficient vanishes. Accordingly, the frequency equation in the case of pervious surface is 


laij| = 0,4, j = 1,2,3, (10) 


where a, j! 8 are same as Ae y8 with r is replaced by a in the eq. (7). In the case of an impervious 
surface the frequency equation becomes 


[bij] =0,i, j = 1,2,3, (11) 
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here, 
b31 = (Q + Rd?) p(p? — k*) Ky (pa), 


b32 = (Q + Rd3)q(q? — k?) Ki (qa), (12) 
b33 = 0, 
where bj; = ai; 1,9 =1,2,7 = 1,2,38. 


4. Numerical results 

For the sake of numerical work, the dissipative coefficient b is taken zero and hence we obtained 
only real phase velocity. To analyze the frequency equations of axially symmetric vibrations of 
poroelastic cylindrical panel, it is convenient to introduce the following non dimensional param- 
eters 


= PB _ _ R _wN 
a1 = F: ao = g, a3 = H> a4 = Fs 


z=(2), y=}, Z= (P, (13) 
p = pi +2p12 + p2, H = P +2Q +R, Ve = 4, 


N 
me 


=e =ü 2 
m= oo C= c= 


In the eq.(13), c is the phase velocity, m is the non-dimensional phase velocity, ka is the non- 
dimensional wavenumber. Employing the non-dimensional quantities in the frequency equa- 
tion, we obtain a implicit relation between non-dimensional phase velocity and non-dimensional 
wavenumber. For the numerical process, three types of poroelastic solids are considered and 
then discussed. Of three poroelastic solids, two are sandstone saturated with kerosene and wa- 
ter, respectively [11, 12] and the third one is bony element. The physical parameters of first two 
materials pertaining to eq. (13) are given in the Table 1. Further, the values of bone poroelastic 
parameters A, N,Q, R and its mass coefficients p;; are computed following the paper [13].The 
values of Young’s modulus and Poisson ratio are taken to be 3 x 10° and 0.28, respectively 
as suggested in [13]. Phase velocity is computed using the bisection method implemented in 
MATLAB and the results are depicted in the figuresl-8. In all the cases, curves are periodic 
and coincide in nature. Figures 1-8 show the plots of non-dimensional phase velocity against the 
wavenumber pertaining to different elastic foundations 10, 20, 50, and 100 respectively. From 
figures 1,3,5,7, it is observed that phase velocity of Material-I values are greater than phase 
velocity of Material-II for both pervious and impervious surface. From figure 1, 3, it is clear 
that phase velocity of Material-I pervious values are greater than phase velocity of Material-I 
impervious surface. Also, it is observed that the phase velocity for both pervious and impervious 
surface of Material-IT is almost steady and constant beyond the wavenumber 1. From figures 5 
and 7 it is observed that the phase velocities of Material-IT both pervious and impervious surface 
is almost constant beyond the wavenumber 5. Figure-2 show the plot of non-dimensional phase 
velocity against non-dimensional wavenumber in the case of bone. From this figure, it is clear 
that phase velocity is same for both pervious and impervious surface beyond the wavenumber 
7. From figure-4, it is observed that phase velocity is constant for both pervious and impervious 
surface beyond the wavenumber 5. From figure 6 and 8, it is clear that phase velocity is same 
for both pervious and impervious surface beyond the wavenumber 3. 
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Figure 1. Variation of non-dimensional phase Figure 2. Variation of non-dimensional 
velocity with the wavenumber at elastic foundation phase velocity with the wavenumber in the 
(K = 10). case of bone at elastic foundation (K = 10). 


—0O— pervious 


—0— Fer, Mat : 
—o— Impervious 


—o— Per Mat II 
—— Imp, Mat! 
—o— Imp, MatII 


aon 
no 


oO 
a 
ON fOAMON 


Phase Velocitty 
Phase V elocity 


eB we Nat shee yay oy 


wavenum ber 


Figure 3. Variation of non-dimensional phase Figure 4. Variation of non-dimensional 
velocity with the wavenumber at elastic foundation phase velocity with the wavenumber in the 
(K = 20). case of bone at elastic foundation (K = 20). 
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Figure 5. Variation of non-dimensional Figure 6. Variation of non-dimensional 
phase velocity with the wavenumber at phase velocity with the wavenumber at 
elastic foundation (K = 50). elastic foundation (K = 100). 
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Figure 8. Variation of non-dimensional 
phase velocity with the wavenumber in the 
case of bone at elastic foundation (K = 


Figure 7. Variation of non-dimensional phase 
velocity with the wavenumber at elastic foundation 


(K = 100). 100). 
Table 1. Material parameters 
Material parameters Materia-I Material-II 
ay 0.843 0.96 
ag 0.065 0.006 
a3 0.027 0.0289 
a4 0.234 0.412 
dy 0.901 0.876 
d2 -0.001 0 
d3 0.1 0.124 
£ 4.869 4.2977 
y 0.998 0.912 
A 3.85 2.126 
5.Conclusion 


Axially symmetric vibrations in poroelastic solid cylindrical panel resting on elastic foundation 
are investigated in the framework of Biot’s theory. Non-dimensional phase velocity against non- 
dimensional wavenumber is computed for three types of poroelastic solids for different elastic 
foundations. Of three poroelastic solids, two are sandstone and third one is bony material. From 
the results, we can infer nature of surface has little influence over the values in the case of bone 
unlike sandstone cylinders. 
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